The characterizing variable for critical point in momentum space 
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The possible experimentally observable signal in momentum space for the critical point, which is 
free from the contamination of statistical fluctuations, is discussed. It is shown that the higher order 
scaled moment of transverse momentum can serve as an appropriate signal for this point, provided 
the transverse momentum distribution has a sudden change when energy increases passing through 
the critical point. A 2-D percolation model with a linear temperature gradient is constructed to 
check this suggestion. A sudden change of third order scaled moment of transverse momentum is 
successfully observed. 
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I. INTRODUCTION 

Heavy-ion collision experiments at the relativistic 
heavy ion collider BNL-RHIC have found evidences that 
a phase transition from hadron gas to quark gluon plasma 
(QGP) has occurred Theoretical studies predict the 
existence of a critical point, or critical end point (CEP) 
— a point which separates the first order phase transition 
at high baryon chemical potential, low temperature from 
the smooth crossover at low baryon chemical potential, 
high temperature 2]. Theoretical estimation strongly 
indicates that, if the critical end point does exist, it is 
within the region of the phase diagram probed by the 
heavy ion collision experiments. Most probably it will 
appear on the quark-gluon phase transition boundary in 
the range of baryon chemical potential of 100-500 MeV, 
which corresponds to heavy ion collisions at ^sjvAf= 5 - 
50 GeV 3]. The two collaborations in RHIC ~ STAR 
and PHENIX — have already begun their procedure of 
low energy scan to observe the critical point 

In order to locate the critical point we must have some 
variable as appropriate signal. This variable should be 
determined by the thermodynamic property of critical 
point. Two distinct properties of the critical point are: 

A, Some quantities have large fluctuations at the crit- 
ical point. A well-known example is the critical opales- 
cence, which appears at the critical point of liquid-gas 
system, indicating the existence of large density fluctua- 
tions. 

B, Some quantities change from zero to non-zero at the 
critical point. For example, in a ferro-magnetic system 
the global magnetization changes from zero to non-zero 
at the critical point. 

These two properties are closely related in thermody- 
namics, but phenomenologically they have different rep- 
resentations and lead to different experimental observ- 
ables. 

In the study of relativistic heavy ion collision the 
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above-mentioned property A has been widely utilized by 
many authors to propose signals for critical point @. It 
is expected that large fluctuations exist at the critical 
point for many variables, such as transverse momentum, 
K/tt ratio, etc. What is important is that these variables 
have to be measured event by event. Due to the limited 
number of particles in a single event, they have the dis- 
advantage of being inevitably contaminated by the sta- 
tistical fluctuations 0, 0, @] ■ There is no reliable method 
to eliminate the statistical fluctuations and it is unclear 
whether any fluctuation signal for critical point can sur- 
vive after eliminating the statistical fluctuations. 

In the present letter we will discuss the possibility of 
finding some signal for the critical point, basing on its 
property B. Such a variable, if exists, will be measured in 
the whole event sample instead of event by event. There- 
fore, it is free from the troublesome problem of statistical 
fluctuations. The statistical errors can be reduced by en- 
larging the statistics, i.e. increasing the number of events 
in the sample. 

The layout of the paper is as follow: In Sect. 2, we 
discuss the moments of final particle momentum distri- 
bution, which may be candidates for such characteriz- 
ing variable. Then, in Sect. 3, a 2-dimensional site- 
percolation model with temperature gradient, which can 
be used to exhibit this characterizing variable, is pre- 
sented and discussed. Finally. The summary and con- 
clusion is presented in Sect. 4. 



II. MOMENTS OF FINAL PARTICLE 
MOMENTUM DISTRIBUTION 

In trying to find a measure basing on the property B of 
critical point it should be noticed that the only quantities 
that are observable in relativistic heavy ion experiments 
are the momentum, mass, charge, etc. of final state par- 
ticles. On the contrary the variables that changes from 
zero to non-zero at the critical point is usually in the 
coordinate space. The aim of this letter is to find out 
a variable in momentum space that can characterize the 
critical behavior in coordinate space. 

Let us ask a question: what information can be ex- 
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FIG. 1: The average transverse-momentum of cj> in A+A colli- 
sions at different colliding energies. Data taken from Ref. 9]. 



tracted from the transverse momentum of final-state par- 
ticles. The answer is: all the available information about 
transverse momentum is contained in its probabiUty dis- 
tribution, or equivalently, in the moments of the distri- 
bution of all orders. Among them are: 

First order moment (pt) gives the average transverse 
momentum, where (•) means average over the whole event 
sample. 

Second order moment (pj) describes the width of the 
distribution. Measuring in each event and combining 
with the first order moment, it gives the event-by-event 
pt-fluctuation cr^ = pi — pf' ^ where and pi are the 
average of pi and pt in each event, respectively. 

Third order moment {p^) indicates the peak of the dis- 
tribution. 

Fourth and other higher order moments describe the 
detailed shape of the distribution. 

As discussed above the event-by-event fluctuation cr^ = 
p1 —pt^ has been suggested by many authors as a possible 
signal for critical point. It has the disadvantage of be- 
ing contaminated by statistical fluctuations. So we turn 
to the flrst order and the higher than second order mo- 
ments. The present available experimental data on (pt) 
of (f) particle show a clear rise and subsequent saturation 
as the increasing of colliding energy, as depicted in Fig.[T] 
This monotonic behavior, if exists also for other identi- 
fied particles, e.g. pions, is encouraging, it indicates that 
possibly an abrupt jump has occurred at a certain energy. 

Let us define the scaled n-th order moment and its 
reciprocal as 
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Dn = {CnY 



(1) 



Both of them are dimensionless and are, therefore, more 
suitable to serve as a characterizing variable. The higher 
than second order scaled moments portray the shape of 
transverse momentum distribution. If when the energy 
increases passing through the critical point there were 
a sudden change in shape of the transverse momentum 
distribution, then the higher order scaled moments of pt 



might have a visible jump. This phenomenon if exists 
can serve as a good candidate of the critical-point signal. 

Here comes another advantage of this approach in com- 
parison with the usual one based on the large fluctuations 
at the critical point. As we know, the flrst round of en- 
ergy scan will be performed with large steps between dif- 
ferent colliding energies. If it is not by occasion that some 
energy used in the first round just locates at the vicinity 
of the fiuctuation peak, we will see nothing in the first 
round and the subsequent scan has to be carried out in 
finer steps over the whole energy range. On the contrary, 
if the higher order scaled pt moments have a sudden rise 
when passing through the critical point, then already in 
the first round of energy scan we will observe a rise of 
these moments, and most probably the critical point is 
located in the region of the moment-rising. Then the 
further scan could be concentrated in this region and the 
critical point, signaled by the abrupt jump of the scaled 
moments, can be catched easily in this way. 



III. 2-DIMENSIONAL SITE-PERCOLATION 
MODEL WITH TEMPERATURE GRADIENT 

Let us use a simple well-defined model that has known 
critical behavior to exhibit the above argument. The 
model is a 2-dimensional site-percolation model[10] with 
temperature gradient. We take a big circle of radius R 
and randomly locate in it n little circles of radius r (r <^ 
R). These little circles will be referred to as cells. If the 
number of cells is large enough, two or more cells may 
overlap. In this case, a connection is built among these 
cells and a cluster is formed. Clusters can be of various 
size. If there is a cluster extending from one side of the 
big circle to the other side, then it is referred to as an 
infinite cluster. The appearance of an infinite cluster is 
considered as the appearance of a new phase. 

The size of the system is deflned as the radius of the 
big circle with that of the small circle as unit: 



r 



(2) 



In the following calculation, we will flx the size of the 
system and assign L the value 500. The model parameter 
is chosen to be the sum of the area of all the small circles 
and the area of the big circle: 



rnir n 



(3) 



In the real physical system the parameter rj corresponds 
to the energy density of the system. 

In Fig. [2] is shown the dependence of the probability 
Poo of the appearance of inflnite cluster on the parameter 
77. It can be seen from the figure that, when 77 is small, 
inflnite cluster almost does not appear, but as rj increases 
when T] goes near to a certain value 77c > the probability 
Poo of infinite cluster increases gradually to unity. The 
bigger the system, the sharper the rising of Poo- When 
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FIG. 2: The dependence of Poo on t). 



L oo, the value of Poo jumps suddenly from to 1. 
The value of r] where this jump appears marks the critical 
point and is denoted by r^c- In our system, r]c is about 
1.12. It is noticeable that all the values of Poo for various 
system sizes intersect at rj — rjc- This is a basic property 
of the percolation model and can be utilized to locate the 
critical value rjc- 

The step- function singularity of Poo [HI is the critical 
phenomenon in the model. It is evidently in the coordi- 
nate space. In order to see how this phenomenon repre- 
sents itself in the momentum space, we add a tempera- 
ture gradient into the model, assuming that the center 
of the big circle is at a higher temperature T^ax — 400 
MeV, and the fringe gets a lower temperature Tmin = 150 
MeV. A linear temperature gradient is assumed along the 
radius of the big circle. 

We further assume that each cluster has arrived at 
thermal equilibrium and all the cells in it have the same 
temperature. The temperature of the ith cluster is de- 
termined by the distance to the center of the big circle 
of its center of mass defined as 



(4) 



where is the distance of the jth cell in cluster i to the 
center of the big circle, and rii is the number of cells in 
cluster I, Thus, the temperature of all the cells in cluster 
i is: 



(5) 



In Fig. [3] are shown the distributions of the number of 
cell at different temperatures. 

From the figures we can see that, as the increase of 77, 
when 77 1.12 a peak appears suddenly in the distribu- 
tion around T ~ 233 MeV, and the height of the peak 
increases very fast. 

The transverse momentum of the cells in a cluster can 
be calculated from the thermal equilibrium momentum 
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FIG. 3: The distribution of cell temperature for different r\. 
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FIG. 4: (Color on line) The distribution of cell transverse- 
momentum for different i). 



distribution, i.e. the Boltzmann distribution 
dA^ N 



2nptdpt 2nmkT 



exp 



2mkT 



(6) 



where k is the Boltzmann constant, and m is the mass of 
the cell. In our calculation, we assume the mass of cell to 
be the mass of pion, equal to 135 MeV. Calculating the 
transverse momentum of each cell from this distribution, 
the pt distributions are obtained and shown in FigH) 

It can be seen from Figdithat the pt distributions for 
rj < 1.15 and 77 > 1.2 fall in two groups with a gap in 
between. This indicates that there is an abrupt change 
in the shape of pt distribution when 77 increases passing 
through the critical point. 

In order to find a variable that characterizes this 
abrupt change of pt distribution, we calculate the mo- 
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FIG. 6: The dependence of D3 on 77. 



ments of pt distribution of first, second and third orders, 
as shown in Fig.[5l It can be seen from the figures that 
both the 1st and the 3rd order moments {pt) and {p^) 
change abruptly at the critical point, while the 2nd or- 
der moment (p^) does not. 

Although the changes of the first and third order mo- 
ments (pt) and (p^) are in the reverse direction, they 
occur at one and the same point, i.e. the critical point. 
The dependence of D3 defined in Eq. (1) on 77, shown in 



Fig. 6, represents a clear step-function shape. In this fig- 
ure the abscissa 77 can be changed to the energy density 
£ and the latter in turn is related to the colliding energy 
Y^snn in heavy ion collision experiments. We see that a 
sudden jump of -D3 as the increasing of ^snn can be used 
as an appropriate signal of critical point in momentum 
space. 

As stated above a strict step function singularity exists 
only in an infinite system. In the real case of heavy ion 
collisions the system is of finite size. The rise of the 
characteristic variable, D^, may not be very sharp and 
may be hard to be recognized. To solve this problem 
we can run the experiments with nuclei of various sizes. 
The intersection of the experimental results will mark the 
place of critical point, cf. Fig. 2. 



IV. SUMMARY AND CONCLUSION 

In this letter we try to utilize the sudden change of 
some variable(s), instead of the large fluctuations, at the 
critical point to find an appropriate variable for locating 
this point. The advantage of such an approach is that it 
is free from the trouble on the elimination of statistical 
fluctuations and that the possible region where the crit- 
ical point locates can easily be fixed already in the first 
rounds of energy scan. Using a site-percolation model 
with temperature gradient we found that the third order 
scaled moment of pt distribution (or its reciprocal) can 
possibly serve for this purpose. 

The model used in this letter is far from a realis- 
tic model of relativistic heavy ion collisions, but the fi- 
nal result is largely model- independent. Checking the 
whole derivation carefully we can see that the essential 
point for to serve as an appropriate signal for crit- 
ical point is the fact that the shape of pt distribution 
changes abruptly at the critical point. The higher order 
scaled moments (or their reciprocal) represent the shape 
of Pt distribution and, therefore, can serve as an appro- 
priate signal. An abrupt change in pt-distribution shape 
when temperature increases passing through the critical 
point is probable for the system produced in relativistic 
heavy ion collisions. The higher order scaled moment or 
its reciprocal is, therefore, a possible candidate for the 
critical-point signal in these collisions. 
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